The Classical Approximation for Real-Time Scalar Field Theory at Finite 
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The use of classical thermal field to approximate real-time quantum thermal field theory is 
discussed. For a Xcj> 4 theory, it is shown that the classical Rayleigh- Jeans divergence can be canceled 
with the appropriate counterterms, and a comparison is made between the classical and quantum 
perturbative expansion, ft is explained why Hard Thermal Loops prevent the same method to work 
for gauge theories. 



I. INTRODUCTION 

Real-time dynamics of quantum fields plays an im- 
portant role in the early universe (baryogenesis, infla- 
tion) and in heavy ion collisions. Some quantities can- 
not be determined within perturbation theory, even af- 
ter possible resummations. An example is the rate of 
sphaleron transitions at high temperature in (extensions 
of) the Standard Model, relevant for baryogenesis. A 
non-perturbative real-time calculation on a lattice, us- 
ing Monte Carlo methods, is complicated, because of the 
complex Boltzmann weight. 

A relatively simple way to do the dynamics is to use the 
classical equations of motion. Depending on the choice 
of initial conditions, this describes an equilibrium or a 
non-equilibrium situation. In this talk I discuss classical 
A</) 4 theory, in thermal equilibrium. This talk is based 
on work done in collaboration with Jan Smit (sec also 
|^| ). More recent work is presented at this Workshop at 
the poster session S. 



This typically leads to a transport or kinetic theory de- 
sciption. As an explicit example, consider real-time Xcj) 4 
theory at finite temperature The low- momentum 

or soft modes are considered to be the degrees of free- 
dom that are of interest, and the high-momentum or 
hard modes as those that can be integrated out. In gen- 
eral, the resulting semi-classical dynamics for the soft 
modes is very complicated. They are subject to noise and 
dissipation, due to the coarse-grained interaction with 
hard modes, and the effective equations of motion are 
non-local in space and time. Especially this non-locality 
makes it difficult to use the resulting equations directly 
for a numerical treatment. 

Therefore I discuss here another possibility to approx- 
imate the dynamics, originally proposed in ||, which is 
related to dimensional reduction (DR) and 3-d effective 
field theories for static quantities: classical thermal field 
theory. 



III. DYNAMICS AND CLASSICAL THERMAL 
FIELD THEORY 
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II. EFFECTIVE FIELD THEORIES AND 
COARSE-GRAINING 

It often happens that a Quantum Field Theory (QFT) 
is too complicated to be solved completely. If there are 
several scales present in the model, a useful approach is to 
construct an effective field theory for the degrees of free- 
dom that are important at the scale under consideration. 
The other 'less important' degrees are integrated out, one 
way or the other. This very general idea has been applied 
successfully in many different physical situations, under 
names as coarse-graining, Wilson renormalization group, 
dimensional reduction, and so on. 

These ideas have also been applied when considering 
the real-time dynamics in a QFT at finite temperature. 



As already mentioned in the Introduction, a relatively 
simple way to deal with the dynamics, is to use the clas- 
sical equations of motions. One simply solves the equa- 
tions, with given initial conditions, and then averages 
over the initial conditions with the Boltzmann weight 
exp — /3H as distribution function. If the hamiltonian in 
the Boltzmann weight is the same as the one that deter- 
mines the equations of motion, this leads to an equilib- 
rium description. For e.g. the classical 2-point function, 
this procedure means the following (x = (x, t)) 



S(x-x') = (0(cc)<Kx')>ci 



(1) 



= z: 



VitV^e-^^ (j){x)(j){x'), 



with the classical partition function 
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In (|l|), is the solution of the classical equations of 
motion 4>{x) = {<p(x), H}, tt(x) = {ir(x), H}, with the 
initial conditions 0(x, to) = </>(x), 7r(x, to) = tt(x). The 
integration over phase space is over the initial conditions 
at t = to, weighted with the Boltzmann weight. 

It is necessary to make two remarks at this point. First 
of all, classical thermal field theory contains the well- 
known Rayleigh-Jeans divergence. Hence, everything is 
formulated with a (lattice) cutoff, to regularize this di- 
vergence. We will show that this cutoff can be removed 
in the end, if the parameters in the classical theory are 
chosen in the correct way. And secondly, when restricting 
to time-independent correlation functions, the canonical 
momenta can be integrated out, and the resulting par- 
tition function has precisely the form of that of a 3-d, 
superrenormalizable field theory, 



z 3-d = I V(f>e 



which has been studied in great detail in DR ||. I will 
come back to DR results at the appropriate places. 

In the remainder of this talk, I will discuss what role 
the cutoff and the Rayleigh-Jeans divergence play and 
what the relation with the dynamics in the quantum the- 
ory is. This will be done in perturbation theory. 



IV. PERTURBATION THEORY IN QUANTUM 
AND CLASSICAL THERMAL FIELD THEORY 



Perturbation theory in the classical theory is obtained 
by writing the field as 4>{x) = 4>o(x) + X(j)i(x) + . . .. Solv- 
ing the equations of motion, order by order in A, gives 

f d 3 k 

<h( x ) = / 77r-vT elk ' X [^( k ) costJ k(i-io) 



(2tt) 3 

7r(k) . , 
H smuj k (t - t )\, 

<h(x) = -\J d 4 x' G*(x - x')^l(x'), etc, 

with u>\ = k 2 + tv? . Here we introduced the retarded 
Green function 

Gl(x - x') = -6{t - t'){{4>{x), <t>{x')}) cX . (2) 

In the unperturbed case, it reads (after spatial Fourier 
transformation) 



G$(k,t-t') =9(t-t') 



sin uj k (t — t') 



(3) 



If one now calculates ([[]) , products of 4>o have to be aver- 
aged with the Boltmann weight, and for the unperturbed 
case this gives 



Sb(M) = J d 



3 xe~ l]i 



*( Mx ) M 0)) cl = T C -^. (4) 



As is well-known, there is a zoo of methods available to 
formulate perturbation theory in real-time quantum ther- 
mal field theory 0. When using the Schwinger-Keldysh 
contour, the quantum field <fi is denoted with <p + , <f>~ , if 
it lives on resp. the upper and the lower branch. A ver- 
sion that is particulary convenient here, is the 'center-of- 
mass/relative' coordinates version, where the basic field 
variables are taken as 



in terms of </> + , <p~ . The reason is that the matrix prop- 
agator in this basis compares directly with the 2-point 
functions that we found in classical perturbation theory. 
Namely, the 2x2 propagator is given by 



iF{x-x') G* m (x-x') 



G(x x ) — ( n A K („ _ J/\ qm q 



with in terms of the original field 

F(x - x') = i(<Mx)<Mx') + <i>{x')ct>{x)), 
G^(x - x') = G^(x' -x) = iO(t - t')([cj>(x), <j>{x>)]). 

Indeed, in the naive classical limit, i.e. fields commute 
and commutators go to Poisson brackets, these become 
(0) and @). 

In the unperturbed case, the relation becomes very ex- 
plicit. The free Gq is given by the classical expression 
(0), and the free Fq is given by 



F (k,t) = [n(uj k ) + -]— — = T 2_^ 



2 Wk 



W„ + Mir 



where n(w) = (e 13 ^ — 1) _1 is the Bose distribution. We 
have written the real-time 2-point function that contains 
all the temperature dependence as a sum over Matsubara 
frequencies Lu n — 27mT, familiar from the imaginary-time 
formalism. Comparing this result with (|j) indicates that 
the classical theory is an effective theory for the n = 
term, just as in DR for time-independent quantities. 

An important ingredient in the setup of perturbation 
theory in the quantum case is the KMS condition || . Us- 
ing the KMS condition to determine the 2-point function 
in temporal momentum space, ensures that the so-called 
vertical part of the contour is taken into account, or in 
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more physical terms, that the system is in thermal equi- 
librium. It turns out that it is possible to derive a KMS 
condition also in the classical theory M, it reads 



-{<t>{x)<l>{x')) cl = T{{ct>{x),<l>{x')}) ch 

It can be used to express S in terms of G R ,G A . Ex- 
plicitly we find in resp. the quantum and the classical 
case the following relation in temporal momentum space 
(k = (k°,k)) 

F (k) = -i(n(k°) + i) (G R (k) - G A (k)) , 
T 



S (k) = -i^(G R (k)-G A (k)) 



Using this, we can work in temporal momentum space 
throughout, which is convenient from a technical point 
of view. 

All the found relations between the quantum and the 
classical expressions are nice, but rather academic until 
we start to calculate perturbative corrections. Hence we 
calculated the 2-point function to two loops and the 4- 
point function to one loop. Both in the quantum as in 
the classical theory, it is straightforward to recognize dia- 
grams and identify 1PI parts. It turns out that there is a 
direct correspondence with the diagrams in the quantum 
theory if one uses the </>i,</>2 basis. I will show this for 
the one loop and two loop setting sun retarded self-energy 
Sh(p). 



One loop self-energy 



In the quantum theory the one loop retarded self- 
energy is given by (see Fig. 1) 
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(5) 



In the classical theory, we find (with a momentum cut- 
off A) 



1 



J R,c\ 



= -A 



d 4 fc AAT 



XmT 
8tt 



(0) 



Compared with the quantum expression, Fq is replaced 
by Sq. The ft term in (^|) is correctly reproduced by 
the classical theory, and the higher order terms in Ti are 
absent, as expected. The term turns up in the classi- 
cal theory as a linear divergence, and not as the thermal 
mass. However, if we use our knowledge from DR, we 
can simply cancel the divergence and put in the correct 
thermal mass by the appropriate choice of classical mass 
parameter, which is written as m 2 = m 2 h — Sm 2 . In DR, 
this is called matching. Note that since this one loop 
diagram is momentum independent, the analysis is not 
very complicated. 



Two loop setting sun diagrams 

The two loop setting sun contribution to the self- 
energy is momentum-dependent and gives in the quan- 
tum theory rise to e.g. Landau damping, because of an 
imaginary part. Hence it is worthwhile to analyse this di- 
agram also in the classical theory. In the quantum theory 
we have two diagrams (see Fig. 2) 
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F (k 1 )F (k 2 )G R (k 3 ) 



-^G fl (fc 1 )G?(fc 2 )G fl (fc 3 ) 

with &3 = p — k\ — k 2 . Note that the second diagram does 
not contain (explicit) temperature dependence. A closer 
look shows that it is indeed subdominant at high T and 
small A. 




FIG. 2. Retarded self energy: two loop setting sun dia- 
grams. The second diagram is absent in the classical theory. 



where I indicated explicitly the Ti dependence. The lead- 
ing term is the so-called thermal mass, w 2 h , which is the 
only Hard Thermal Loop (HTL) contribution in scalar 
field theory. It needs to be resummed in order to have a 
true perturbative expansion (in A 1 / 2 ). 




FIG. 1. Retarded self energy: one loop diagram. The full 
line represents Fq resp. So and the dashed- full line Gq. 



It turns out that in the classical theory we only find 
the dominant contribution, i.e. the classical counterpart 
(Fq — > So) of the first diagram 



£fl UI ci(p) 



iA 2 
2 

A 2 T 2 



(i 



d A k Y d 4 k 2 
(2tt) 4 (2tt) 4 
d 3 h d 3 k 2 
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' T 



d£l w(~p, fi) 



2tt p° + it + n ' 
where we introduced the classical scattering integral 



(7) 
(8) 
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»<P,n) = ir£/n 



{±} 3=1 
4 



T 



(27r) 3 2w kj . Wk 3 



x(27r) 4 <5(fi±w kl ± w k2 ±Wk 3 ). 

The sum is over all +'s and — 's. The first term ((?]) is in- 
dependent ofp°, real, and logarithmically divergent. It is 
actually the result obtained in DR. Hence the divergence 
is canceled again by the appropriate choice of classical 
mass parameter, just as in the one loop diagram. The 
divergences we have encountered by now are the only 
ones in DR, because of superrenormalizability. A closer 
look at other diagrams in the classical theory makes it 
plausible that these divergences are also the only ones 
in the general time-dependent classical case. Hence the 
classical divergences are completely under control. The 
second term (g) is finite. It contains all the p° i.e. time 
dependence, and it has an imaginary part. 

It is instructive to compare the imaginary parts in the 
quantum and classical case. Using the notation E9] (for 
p Q > 0) 



Im £™(p) 



-9i(p) ~ 92(p), 



for both the quantum and classical case, we find, in the 
quantum theory 



o /T f d 3 kx d 3 k 2 ra kl 



A^ 

n k2 77 k3 „ c/ o \ 
x 2nd{p - uj kl - w k2 - w k3 ), 

L0k 2 Wk 3 

32 7 (27r)' ;i (2-7r) d w kl 

77 k2 ^k 3 r» c/" fl , \ 
x 27rd(p + w kl - ^k 2 - Wk 3 ), 

w k 2 w k 3 

and in the classical theory 

AV f d 3 ^ d 3 k 2 T 3 
9i AP) -wY] (2n) 3 ( 27 r) 3 «< 

x2ir5(p t) - w kl - w k2 - ), 
A 2 p° f d 3 k! d 3 k 2 T 3 



A p' r d 3 ki d 6 h 

92,d(p) - J (2V) 3 7>o 



O 3 

x2ir5(p° + uj kl - uj k2 - w k3 ), 

5i,(ci)(p) represents three body decay, and <72,(ci)(p) Lan- 
dau damping. The classical expressions are indeed the 
leading order expressions from the quantum theory. The 
reason is that, after replacement of the Bose distributions 
by the classical distributions, the resulting integrals are 
still finite. 

To demonstrate this very explicitly, consider the on- 
shell plasmon damping rate (Ojjlj]. It is given by 



7 = 



-ImSr(0,m) _ 32(0, m) 



It turns out that, to leading order in A and T, 52(0,771) 
equals 32.ci(0, m), and we find (to leading order) 



7ci = 7 = 



A 2 r 2 

15367T777 



Using that m 777th, because of resummation in the 
quantum theory and because of matching in the classi- 
cal theory, gives for the plasmon damping rate (again to 
leading order) 



7ci = 7 



\VXhT 
128%/6tt' 



This example shows in a very explicit way that the clas- 
sical theory, with the appropriate choice of parameters, 
indeed approximates the quantum theory. 

Let me end the discussion on scalar field theory with 
the remark that also for the 4-point function the leading 
order expressions are reproduced by the classical theory. 



V. CLASSICAL THERMAL GAUGE THEORY 
AND HARD THERMAL LOOPS 

Of course, one of the reasons the scalar theory was 
studied, is because of possible applications to the dynam- 
ics of hot gauge theories. We know that the DR approach 
has been succesfully applied to study e.g. the electroweak 
phase transition, and the presence of (non-)abelian gauge 
fields did not pose any fundamental problems. 

It turns out that, concerning the dynamics, there are 
such problems 1 12 lj| . The reason for this are the Hard 
Thermal Loops in gauge theories E3). In contrast to 



the scalar theory, HTL's in gauge theories are momen- 
tum dependent, in a complicated non- analytical way. As 
a typical example, consider the longitudinal part of the 
gauge boson self-energy in e.g. scalar electrodynamics. 
In the HTL approximation, it is given by 

\ p V 2 p p - pi 



with the prefactor 



77/ 



7T 2 , 

7T ./o 



dkkn(k) 



e 2 T 2 
6h 



where I indicated again the h dependence. This term 
gives rise to e.g. the Debye screening mass, m 2 D — 
F(0,p -> 0) = 27T7 2 = e 2 T 2 /3h. 

In the classical theory, the Bose distribution n(u>) is re- 
placed by the classical one T/ui, which leads to a linearly 
divergent prefactor 



= — - / dkk 1 - = 



3 2 TA 



2m 



2m 



where we used a momentum cutoff for simplicity. This 
is similar to what happens in the scalar theory (compare 
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(0) and (||)). However, now it is not possible to cancel 
the divergence with a local mass counterterm, because of 
the complicated momentum dependence. In the DR ap- 
proach for static quantities, these problems do not arise, 
because for p° — 0, -F(0,p) = 2m 2 , i.e. momentum inde- 
pendent. Hence the linear divergent Debye mass in the 
effective 3-d theory can, in DR, simply be canceled by a 
mass counterterm. 

We can conclude that HTL's in gauge theories make 
a straightforward use of classical thermal gauge theory 
to approximate the dynamics questionable. In the sev- 
eral proposals that exist to incorporate HTL effects in 
a classical-like theory, new local degrees of freedom are 
added to make up for them in one way or the other 
|^2| , ^5| -[l7| . The idea to add classical particles (instead 
of fields) has been succesfully implemented numerically 



VI. CONCLUSIONS 

To summarize, we have shown how classical thermal 
field theory can be used to approximate real-time quan- 
tum field at finite temperature, for the A(/f> 4 case. Instead 
of explicitly integrating out hard modes to construct an 
effective theory for the soft modes, we showed that it is 
possible to take all the modes into account. The resulting 
Rayleigh-Jeans divergence can be dealt with in a straight- 
forward manner, namely by using countcrterms that are 
dictated by dimensional reduction. Furthermore, by us- 
ing real-time perturbation theory, we have shown that the 
classical theory approximates the quantum one if also the 
finite part of the classical parameters are chosen accord- 
ing to the DR matching rules. 

Essential is that the HTL effects in the quantum theory 
can be easily incorporated in the classical theory: only 
the thermal mass has to be put in. This is also the reason 
why the same prescription does not work for gauge the- 
ories. Here HTL's are momentum dependent and in the 
classical theory, they give rise to divergences that cannot 
be canceled with local counterterms (i.e. in 3 + 1 dimen- 
sions; the classical approximation does work for gauge 
theories in 1 + 1 dimensions, because here the classical 
contribution is dominant and finite, and the HTL con- 
tribution is subdominant |lS||). Adding new degrees of 



freedom to represent the HTL contributions seems to be 
a possible way out. 
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